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On the linear stability of collisionless microtearing modes 
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Microtearing modes are an important drive of turbulent heat transport in present-day fusion plasmas. We 
investigate their linear stability under very-low collisionality regimes, expected for the next generations of de- 
vices, using gyrokinctic and drift-kinetic approaches. At odds with current opinion, we show that collisionless 
microtearing instabilities may occur in certain experimental conditions, particularly relevant for such devices 
as reversed field pinches and spherical tokamaks. 
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Microtearing modes (MTMs) are short-wavelength 
electromagnetic instabilities driven by electron temper- 
ature gradients, investigated since the 1970s for their 
relevance to heat transport in fusion plasmas^. Inter- 
est toward them has strongly revived in the last decade 
thanks to the widespread adoption of sophisticated gy- 
rokinctic codes running on powerful computers, allow- 
ing the study of MTMs in realistic geometries, including 
tokamaks, spherical tokamaks, and reversed field pinches 
(RFPs)2r— . Recently the first nonlinear simulations of 
MTMs^ were able to show that these modes can be re- 
sponsible for a large fraction of turbulent heat transport 
in tokamaks. Apart from the plasma core region, MTMs 
are thought to be relevant also in the edge as players 
regulating the heat transport, and possibly the pedestal 
evolution between ELMs^. 

Even more important is to assess how MTMs can 
impact future plasma conditions, in the quest towards 
reactor-grade scenarios. To a large extent, this amounts 
to investigating hotter and thus less collisional plasmas. 
The reference works addressing the role of collisionality 
v are Drake et a&, and Gladd et a&. Their conclusion is 
that - all the other parameters fixed - the growth rate 
peaks at a finite value of v and decreases down to nega- 
tive values as v — > 0. 

This basic picture held until nowadays, just marginally 
modified by subsequent studies by Catto and Rosen- 
bluthi^ and Garbet et aiii, and is apparently supported 
by gyrokinetic simulations^. Very recently, on the other 
hand, several papers have been accumulating a consis- 
tent amount of evidence in favor of the opposite view: 
MTMs may actually be destabilized even under vanish- 
ingly small-^ conditions^— . However this disagreement 
is a consequence of the fact that a generic collisionless 
regime does actually not exist: v is not the only impor- 
tant non-ideal mechanism; its role may be played, for 
example, by particle inertia. These mechanisms were 
cither unaccounted for in earlier works, or the plasma 
conditions addressed were not suitable for their develop- 
ment. 

The purpose of this Letter is to investigate the occur- 
rence of collisionless MTMs, considering previously unac- 
counted destabilizing mechanisms/regimes. They include 
the magnetic drifts in the kinetic equations, and the mu- 
tual role of particle density and magnetic shear profiles. 



At first we will make an assessment on MTM linear stabil- 
ity using the gyrokinctic code GS2i£ in the experimental 
scenarios more prone to MTM turbulence, in particu- 
lar for the reversed field pinch geometry. Afterwards we 
will validate such results with a slab drift-kinetic model, 
and justify our conclusions as relevant in particular for 
high magnetic drift configurations. Besides yielding an 
independent cross-check of the results, the drift-kinetic 
model, although less accurate quantitatively, provides a 
deeper insight into the physical mechanisms behind the 
mode dcstabilization. 

Gyrokinetic approach. The gyrokinctic equation is 
solved in a flux-tube domain by means of the electromag- 
netic code GS2 in RFP geometry^. Unlike the tokamak 
configuration, the RFP is characterized by a fast wind- 
ing of the magnetic field in the poloidal direction, which 
causes the safety factor profile q to be <C 1 , and to vanish 
in the very edge of the plasma. This reflects on a differ- 
ent ordering of the drift frequency tud in the gyrokinctic 
equation. Given B the magnetic field strength (normal- 
ized to the average field on the flux surface), with Bg 
and B,j, the poloidal and toroidal component, the gradi- 
ent uvs and curvature uiu terms in u>d oscillate around 
—B' (hereafter ■' = d ■ /dr) and Bj/r respectively as a 
function of the poloidal (ballooning) angle 9, see Ref. HH. 
Thus the average normalized uvb and ujk are not of the 
same order as the inverse aspect ratio a/R, like in the 
tokamak, but are typically larger than 1. 

For our study we start from a reference internal- 
transport-barrier scenario of RFX-mod, already inves- 
tigated in Ref. [3: experimental values of safety factor 
q ~ f/fO, magnetic shear s = rq'/q ~ —0.6, normal- 
ized logarithmic density gradient a/L n = —an'/n — 0.2, 
with a torus minor radius; for the ion/electron temper- 
ature gradients we assume a/L^ = (in order to limit 
the number of potential instabilities, this quantity is ac- 
tually unknown in RFX-mod) and a/Lx e — 4, respec- 
tively. The plasma f$, the parameter ruling the electro- 
magnetic effects in the gyrokinctic equation, is artificially 
increased with respect to the experimental value (/3 e ap- 
proximately 0.7%), so as to work well above the stability 
threshold. Particle collisions are included by means of 
a classical Lorentz operator, fitted to describe pitch an- 
gle scattering and trapped-passing particle interactions; 
the momentum exchange between species (electrons and 
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ions), together with the energy scattering and slowing 
down processes, is not considered in GS2. For the nu- 
merical simulations we include fluctuations in the paral- 
lel vector potential An. The electrostatic potential <j> is 
usually retained, although it may be artificially switched 
off in some cases, so as to assess its impact upon final 
results. Throughout this section, v is normalized to its 
experimental value at the plasma conditions considered, 
h> cxp ~ 0.4wth,i/a, with vth.i ion thermal speed. Con- 
cerning the resolution of our calculations, the longitudi- 
nal grid must be large enough to resolve the elongated 
structure of the MTM eigenfunctions, \9\ < 607T with 32 
grid points per 27r-period; for the velocity space we typ- 
ically set 12 untrapped pitch angles and 16 values of en- 
ergy, having performed convergence tests with increased 
number of grid points on a subset of cases at the lowest 
collisionalities. 

For all such parameters the fastest growing instability 
turns out to be of MTM type, in the range k y pi < 1. 
Sometimes, another branch with opposite parity appears 
at higher wavenumbers, due to the passing electron re- 
sponse, similar to that described in Ref. Il9l . The mi- 
crotearing nature of the modes is clearly evident by look- 
ing at the mode structure (odd symmetric and very elon- 
gated in <fi, even and localized in Aii), at the sign of the 
real frequency (corresponding to a propagation in the 
electron diamagnctic direction), and at the trend of the 
growth rate as a function of electron temperature gradi- 
ent and plasma /3. The growth rate definitely docs not 
vanish when decreasing collisionality, as shown in Fig. [1] 

The most relevant difference between the slab model 
used in Refs. [H and|H and the RFP geometry is the com- 
plexity of the LOd drift in the latter. In the following we 
investigate its role in destabilizing MTMs, having defined 
f'VBifk € [0,1] as parameters artificially weighting the 
fraction of wvb and uik respectively, keeping /vs = fk 
for the moment. In Fig. [T] we show the growth rate of 
the fastest growing MTM for a flat-density (a/L n = 0) 
high-beta plasma. While the curves have a maximum 
at a given is, the growth rate is non- vanishing also for 
v — > 0, except for the case without <f> fluctuations and 
no drifts. In agreement with Ref. |U, retaining the elec- 
trostatic potential is always destabilizing. Large growth 
rates for /vs = /fc = 1 generally correspond to upshiftcd 
"t(k y ) spectra, peaked for k y pi > 1/2; otherwise, the peak 
occurs around k y pi = 0.1 — 0.3. 

It is instructive to quantify how the single drifts sep- 
arately modify the destabilization of MTMs, Fig. [2j For 
this study we assume 4> ^ 0, a higher wavenumbcr 
and plasma beta, (k y pi = 0.5, j3 e = 0.08), and a very 
low collisionality, v = 10~ 6 . The contribution of the 
two drifts turns out to be comparable; they separately 
provide a monotonic increase of the growth rate. For 
fvB = = fk, at this value of k y pi, MTMs are stable. 
Note that this is due to the mentioned downshift of the 
cutoff of the spectrum ~f{k y ) as fvB,fk —> 0. 

Up to this point we may claim that devices featur- 
ing high drifts arc particularly prone to MTMs. Besides 
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FIG. 1. Growth rate as a function of collisionality, with and 
without electrostatic fluctuations, with (/vs = fk = 1) and 
without (/vs = fk = 0) ojvs and u>k drifts. Relevant param- 
eters: a/L T<! = 4, a/L n = 0, s = -0.65, q = 0.12, fl. = 0.05, 
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FIG. 2. Growth rate as a function of fractions of ojvs and ujk 
(a) with some representative eigenfunctions of electrostatic 
potential (b) and parallel vector potential (c), for v = 10~ 6 . 
Relevant parameters: a/LT c = 4, a/L n = 0.2, s = —0.65, 
q = 0.12, /3 e = 0.08, k y pi = 0.5. 



drifts, however, we identify another destabilizing mech- 
anism, which arises due to the mutual balance between 
magnetic shear, temperature and density gradient. In 
Fig. Owe show the growth rate as a function of magnetic 
shear s and logarithmic density gradient a/L n for a low 
collisionality plasma, v = 10~ 3 , and with /vs = fk = 0. 
Furthermore, the trapped particle fraction is set to 0, so 
as to exclude further destabilizing contributions^. We 
find the existence of a maximum of 7 in the a/L n < 
half-plane: positive density gradients are increasingly 
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FIG. 3. Growth rate vs. magnetic shear and density gradient 
for v — 10 -3 (a.u.), including electrostatic fluctuations (a). 
The section § = 1.4 is shown in (b). The same section for 
v = 1 is in (c). Relevant parameters: a/Lr c = 4, q = 0.12, 
p e = 0.05, fvB = = / fe , fc„p< = 0.2. 



stabilizing MTMs for each value of the magnetic shear, 
Fig. [3]-a/b, whereas slightly negative gradients provide 
destabilization. An analogous behaviour of j(a/L n ) has 
been recently discussed in Ref. for (collisional) spher- 
ical tokamak plasmas. On the other hand, for a fixed 
density gradient the growth rate has a maximum in s. 
The appearance of this maximum is not unexpected, 
since a similar dependence on the parameter L n / L Sl with 
L s = qR/s magnetic shear length, was already encoun- 
tered in Ref. and again in Ref. the difference here 
is that we find a positive 7 also for v <C 1. Notice that, 
at this level of collisionality, the only way to destabilize 
MTMs is to include fluctuations. The same study per- 
formed keeping = identically provides MTM stability 
for every couple (a/L n ,s). Increasing collisionality fur- 
ther destabilizes the mode, as expected: in Fig. [3}c the 
function j(a/L n ) is shown for v = 1. While the roles 
of the density gradient and of the magnetic shear are 
confirmed, we notice that the maximum in a/ L n is mov- 
ing to much lower (unphysical) values, and that, for this 
collisionality, the modes can be destabilized even with 
= 0. These results do not contradict Refs. [H and [H, 
since we are in a different regime: those papers addressed 
the case of a plasma featuring strong temperature gra- 
dients and more moderate but still notable density gra- 
dients, 77 = L u /Lt s > 1) £ = L n /L s -C 1, which is a 
fairly reasonable picture of a device centrally heated and 
fuelled. On the other hand, while a reactor-grade device 
must feature relevant temperature gradients, density pro- 
files could be rather flat in the absence of relevant central 
particle sources, as happens in devices without neutral 
beams. 



Drift-kinetic approach. We solve the Ampere equation 
and the quasi-neutrality equation in slab geometry, as 
extensively explained in Refs. H andU: 
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ijjcr — x (p) (1) 
— Xxipa (u>A\\ — x ip) (2) 



where \ = £ 2 (c Q /c) 2 , ip = {S/'i^^uj^/c) 2 , (p = 
4>^{c/ui), and •' = d ■ jdx. The parallel conductivity 
a is defined in the original papers^. Lengths are nor- 
malized to the ion thermal Larmor radius p, and fre- 
quencies to the electron diamagnctic frequency ; x is a 
radial-like coordinate orthogonal to the equilibrium mag- 
netic field Bqz, x — being the resonance position of the 
mode; k y and k\\ = k y x/L s are respectively the compo- 
nent of the wavenumber of the mode orthogonal and par- 
allel to the magnetic field; c,c a ,Ui the light, Alfven and 
ion thermal speed, uj pe the plasma frequency. Equations 
([T]) and @ are to be solved with the boundary condi- 
tions satisfying the symmetry of the MTM: A'(0) = 0, 
0(0) = 0, A'(oo)/A(oo) = 0'(oo)/0(oo) = -ky. This 
yields four real- valued eigenvalue problems for the quan- 
tities u> = L) r + 17 (the eigenfrequency of the mode) and 
Re(0'(O)), Im(0'(O)), that are determined via Newton- 
Raphson method. Our solver of Eqns. (1 1 12[) has been 
first benchmarked against a selected set of results picked 
up from Refs. [H andH. 

According to the previous section, let us first ad- 
dress the issue of the drifts. The slab geometry 
by definition cannot accommodate curvature drifts, 
but V.B ones can be included heuristically: uvs ~ 
— (m e u 2 L /2eB ) {x/L 2 s )y + 0(x 2 )z. At this stage, u VB 
is defined apart from numerical factors of order unity; 
it can be conveniently written in dimensionless form as 
ays = f\/BxL n /L s , where the factor /vb has been de- 
fined above, and is used to modulate the effect of the 
V.B drift. The collisionless linearized drift- kinetic equa- 
tion becomes 



ui — ik\\u\\ — ikyUys] f = 
eE^fo/T - (ikyc/ B )(cf) - A\\U\\/c)f , 



(3) 



where /o is the equilibrium electron distribution func- 
tion, and / its perturbed part. The first moment over 
the velocity yields the perturbed current used in the 
Ampere and quasi-neutrality equations. Indeed, by in- 
creasing /vb the drift term becomes more and more im- 
portant: an approximately linear increase of 7 with /vb 
is featured, as is shown in Fig. |4l The physical picture 
is that non-ambipolar drifts, such as V B and curvature, 
provide the current across the magnetic field needed for 
reconnection. 

We focus then on the role of plasma gradients, with 
no longer consideration of the drifts, /vb =0; here the 
drift-kinetic formalism is expected to be more appropri- 
ate, and allows us to gain insight into the results of Fig. [3] 
Taking simultaneously the limits 77, £ — > 00 and v — > 
considerably simplifies the equations above. Setting for 
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FIG. 4. Growth rate versus VB drift (for the definition of 
fvB see the main text). Other parameters are: v = 0, k y = 
0.05, £ = lO -2 ,?? = !,£ = 0.1 
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FIG. 5. Real (red squares) and imaginary (blue dots) part of 
the eigenfrequency versus collision frequency. Other parame- 
ters are k y = 0.05, = 10~ 2 ,?7 = £ = 6. 



further handiness rj = £, in these limits the conductiv- 
ity a shrinks within a layer of width Ax ~ p e around 
x = 0, a physical consequence of the finite electron in- 
ertia. Being ^(0) = 0, the scalar potential is very small 
for \x\ < p e , hence its correction to the r.h.s of Eq. ^ 
is negligible (at variance with the gyrokinetic approach, 
where <fi plays a more important role). Eqns. ([1]) and ([2]) 
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FIG. 6. Real (red squares) and imaginary (blue dots) part 
of the eigenfrequency versus rj = £. Other parameters are 



v = 10 ,/3 = 10" 
reference. 



, k y = 0.05. 
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decouple, and we can fix u> by just solving the former. 
Inspection of this simplified equation shows that oj and rj 
do not enter individually but through their combination 
ujjrj = A, which becomes now the proper eigenvalue. We 
can now compute A, finding Im(A) > 0, which implies 7 
to be linearly growing with rj for r) 1 (and the same 
dependence holds for ui r ). 

Quantitative results are shown in Fig. [5j where u is 
shown as a function of the collisionality for fixed values 
of the other parameters, and rj = £ ^> 1: a finite 7 > is 
definitely recovered for v — > 0. Fig. [5] shows the depen- 
dence of u on 77 = £ for v -C 1: the linear trend expected 
on the basis of the previous paragraph is found. Hence, 
also in this regime gyrokinetic results are recovered qual- 
itatively. 

Concluding remarks. The linear stability of MTMs 
represents a more multifaceted problem than previously 
thought; a unique "collisionless regime" cannot be sin- 
gled out; rather, parameter regions with different stabil- 
ity properties can be identified depending on mechanisms 
such as particle inertia and non-ambipolar drifts. In par- 
ticular the intrinsic high magnetic drifts of reversed field 
pinches and spherical tokamaks make such devices more 
prone to microtearing turbulence. 
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